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Abstract 



The paper presents a study of kinematic axisymmetric mean-field dynamo 
models for a case of the meridional circulation with a deep-seated stagnation 
point and a strong return flow at the bottom of the convection zone. This kind 
of circulation follows from mean-field models of the angular momentum balance 
in the solar convection zone. We show that it is possible for this types of merid- 
ional circulation to construct kinematic dynamo models that resemble in some 
aspects the sunspot magnetic activity cycle. The dynamo model includes turbu- 
lent sources of the large-scale poloidal magnetic field production, due to kinetic 
helicity and a combined effect due to Coriolis force and the large-scale current. 
In these models the toroidal magnetic field, which is responsible for the sunspot 
production, is concentrated at the bottom of the convection zone, and is trans- 
ported to low-latitude regions by the meridional flow. The meridional component 
of the poloidal field is also concentrated at the bottom of the convection zone 
while the radial component is concentrated in near polar regions. There are some 
issues which, perhaps, are resulted from the given meridional circulation pattern 
and the distribution of the magnetic diffusivity inside convection zone. In partic- 
ular, in the near-equatorial regions the phase relations between the toroidal and 
poloidal components disagree with observations. Also, we show that the period 
of the magnetic cycle may not always monotonically decrease with the increase 
of the meridional flow speed. Thus, for the further progress it is important to 
determine the structure of the meridional circulation, which is one of the critical 
properties, from helioseismology observations. 
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Introduction 



The widely accepted paradigm about the nature of the global magnetic activity of 
the Sun assumes that meridional circulation is an important part of the dynamo processes 



, 1995; 


Durney , 


1995; 


Choudhuri 



& Dikpati, 1999; Miesch et al., 2010). The current flux-transport and mean-field dynamo 



models, (e.g., Dikpati & Charbonneau, 1999; Guerrero & Munoz, 2004; Bonanno et al. 



2002), typically employ an analytical profile of the meridional circulation pattern, which has 
to satisfy a mass conservation equation and the relevant boundary conditions. One of the 
basic features of this profile is that the circulation stagnation point is close to the middle 



of the convection zone, (e.g., Bonanno et al., 2002; Dikpati et al., 2004). However, such 



ad hoc models of the meridional flow have no support from the mean-field theory of the 



angular momentum distribution in the solar convection zone (Kitchatinov & Riidiger, 1999 



Kitchatinov & Olemskoy, 2011), also see, Garaud & Acevedo Arreguin (2009). The theory 



predicts the meridional circulation pattern with nearly equal amplitudes of the flow velocity 
at the bottom and top of the convection zone. The stagnation point of this flow is close to 
the bottom of convection zone, near O.75i? , and the circulation is concentrated near the 



convection zone boundaries. Rempel(2005; 2006) obtained a similar meridional circulation 



profile with a deep stagnation point and used it to construct a nonlinear dynamo model. 

The physical mechanisms of the strong deviation of the meridional circulation pattern 
from the simple analytical models are discussed in the recent paper by |Kitchatinov fc ~ 



Olemskoy (2011). Here, we briefly summarize their main arguments. The distribution of 
large-scale flows in the balk of convection zone is close to the Taylor-Proudman balance. 
However, this balance is violated near the boundaries. This results in a concentration of the 
velocity circulation in the Eckman layers near the bottom and the top of convection zone 



(Durney, 


1999; 


Miesch et al., 


2006; 


Brun et al. , 


2010) 
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Our goal is to investigate how the meridional circulation with a fast return flow at 
the bottom of the convection zone can affect solar dynamo models. We construct a series 
of kinematic mean-field dynamo models that employ the meridional circulation pattern 



suggested by Kitchatinov & Olemskoy (2011). These dynamo models include the turbulent 
generation of the magnetic field due to the kinetic helicity ( a-effect), the combined effect 
of the Coriolis force and large-scale current (fi x J-effect), and the toroidal magnetic field 
generation due to the differential rotation (^-effect). Following Krause and Radler(1980) 
these models can be classified as a 2 8Vt dynamo. Our approach is to investigate conditions 
of the dynamo instability for this type of the meridional circulation and determine the 
basic properties of the dynamo solution at the instability threshold. This is a kinematic 
dynamo problem. The next section describes the formulation of the mean-field dynamo 
model, including the basic assumptions, the reference model of the solar convection zone, 
and input parameters of the large-scale flows. Section 3 presents the results and discussion. 
The main findings are summarized in Section 4. 



2. Basic equations 

2.1. Formulation of model 

The dynamo model is based on the standard mean-field induction equation in perfect 
conductive media (Krause and Radler, 1980): 

— = Vx (£+UxB), (1) 



where S = u x b is the mean electromotive force, with u, b being the turbulent fluctuating 
velocity and magnetic field respectively; U is the mean velocity. General expression for £ 



was computed by Pipin (2008) (hereafter P08). Following Krause and Radler (1980) we write 
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the expression for the mean-electromotive force as follows: 



Si = (aij + 7^) B - r/ijkVjBk 



(2) 



where, tensor a^j represents the turbulent alpha effect, tensor 7^ describes the turbulent 
pumping, and the 77^ term describes the anisotropic diffusion due to the Coriolis force and 



the fix J effect (Radler, 1969). 



We consider a large-scale axisymmetric magnetic field, B = e^B + V x 



where 



B(r,0,t) is the azimuthal component of the magnetic field, A(r,9,t) is proportional to the 
azimuthal component of the vector potential, r is radial coordinate and 6 - polar angle. The 
mean flow is given by velocity vector U = e r U r + eoUo + e^r sin0f2, where Q (r, 6) is the 
angular velocity of the solar differential rotation and f7 r (r, 0), Uq{it,Q), represent velocity 
components of the meridional circulation. The mean-field magnetic field evolution of is 
governed by the dynamo equations, which follow from Eq.Q: 



OA 
~dt 
dB 



. ac , UesmOdA TT OA 

rsmdt^ H U r — 

r oji or 



(3) 



sin^ fdVLdA dttdA\ d{rU r B) + sin 9 dU e B + ldrSg + sm0d£ r 
dt r \dr dji dji dr J dr r dji r dr r dji 

We introduce the free parameter C v to control the turbulent diffusion coefficient (see 
Appendix). A solar- type dynamo model can be constructed not only with the a-effect as 



a prime turbulent source of the poloidal magnetic field generation flStix , 1976; Pipin & 



Seehafer, 2009; Seehafer & Pipin, 2009). The exact mechanism of the large-scale poloidal 



magnetic field production on the Sun is not known. After Parker (1955), it is belived that 



the a-effect (associated with cyclonic convection) is the most important turbulent source 
of the poloidal magnetic field generation on the Sun. The a-effect particularly important 
role for the dynamo because it describes a leading term of the Taylor expansion of the 
turbulent mean-electromotive force in terms of the scale separation parameter. In addition, 
the mean-field theory predicts the magnetic field generation effects due to the interaction 
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of the Coriolis force (fi x J-effect) and the differential rotation (W x J-effect, where W is 



the large-scale velocity shear) with large-scale electric current (see, |Radler , 1969; Krause & 



Radler, 1980; Rogachevskii & Kleeorin, 2003). The role of these mechanisms for the solar 



dynamo is not well understood though it was shown (Pipin & Seehafer, 2008, 2009) that it 



is possible to construct the solar-type dynamo models including these effects. Our model 
includes the a- and (lx J effects. We introduce free parameters C a and Cf 1 ^ to control 
their strength. 



We use the solar convection zone model computed by Stix (2002), in which mixing 



length £ olmlt |A(p>| , where A (p) = Vlog p is the inverse pressure scale height, and 
olmlt = 2. We confine the integration domain between O.712i? and O.972i? in radius. It 
extends from the pole to pole in latitude. The differential rotation profile, Q = f^o/h (x, /x), 
x = r/i? , /i = cos# is a modified version of the analytical approximation to helioseismology 



data proposed by Antia et al. (1998), see Figure la 



The meridional flow is modeled in the form of two stationary circulation cells, one in 
the northern and one in the southern hemisphere, with a poleward motion in the upper 
and equator- ward motion in the lower part of the convection zone. Following |Kitchatinov 



& Olemskoy (2011), the meridional circulation velocity components were approximated via 



the orthogonal Chebyshev polynomial decompositions: 

3 



U e = 3[/ o sin^cos0^4 n) r n (O, 

71=0 

3 

U r = U (3cOS 2 6-l)(l-e)J2 C r n)T n(0 



71=0 



where 



2^y 0C q X 



(5) 
(6) 

(7) 



X e X\) 

Here, X\) e are the radial boundaries of the integration domain. In our case, x^ = 0.712 and 

» » 



x e = 0.972. The coefficients c K e , c K r are given in Table 111 Parameter Uq controls the speed 
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Fig. 1. — Internal parameters of the solar convection zone: a) the contours of the constant 
angular velocity are plotted for the levels (0.75 — 1.05)fio with a step of 0.025f2oj = 
2.86 • 10 -7 5 _1 ; b) the vector field of the meridional circulation with |U| measured in units of 
Uq. c)the meridional component of circulation at 6 = 45°; d) turnover convection time r c , 
the background turbulent diffusivity rj^\ RMS convective velocity U' c \ e) the radial profiles 
of the a-effect components at 9 = 45°; 
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of the meridional circulation. Kitchatinov & Olemskoy (2011) obtain Uq « 16 ms 1 . 



The geometry of the meridional flow is illustrated in Figure [T](b) . Figure [T](c) shows 
the latitudinal component of the circulation in units Uq . 

The boundary conditions represent a perfect conductor at the bottom and the potential 
magnetic field configuration outside the domain. 



2.2. Method of solution 



We investigate the linear stability of the dynamo equations ( 3j4 ) and determine unstable 
dynamo modes. Then we construct linear dynamo solutions using the corresponding eigen- 



functions. Our approach to solve the linear problem was described in details by Pipin 



& Seehafer (2009) and Seehafer & Pipin (2009). We use a Galerkin method, expanding 



the magnetic field in terms of a basis that satisfies the boundary conditions (Boyd, 2001 
Livermore & Jackson, 2005). The system of Eqs. ^ and ([!]) has exponentially growing or 
decaying solutions, which we represent in the form 



A(x,e,t) = e at J2J2 A ^ sines n A nl(OPL(cose), 



(8) 
(9) 



where 5nm (£) and (£) are linear combinations of Legendre polynomials, and is 
the associated Legendre function of degree m and order 1. These expansions ensure the 



n 





1 


2 


3 


>) 


-0.13432(5) 


- 0.40473(6) 


- 0.02170(3) 


- 0.10718(5) 


>) 


- 0.0681469(4) 


- 0.006839(4) 


-0.032516(1) 


- 0.0027(4) 



Table 1: The coefficients for the meridional circulation profile components given by Eqs. g|6[ ). 
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regularity of the solutions at the poles 9 = and 9 = tt. The integrations in radius and 
latitude, which are needed for calculating the expansion coefficients A nm and B nm) were done 
by means of the Gauss-Legendre procedure. The eigenvalue problem for determining the 
eigenvalues, a, and the associated eigenfunctions was solved by using the LAPACK software. 
There are two types of dynamo eigenmodes: 1) modes with a symmetric distribution of 
the toroidal component B and antisymmetric distribution of the poloidal component A, 
relative to the equator, called here "S-modes", and 2) vice-versa antisymmetric al modes, 
called "A-modes". We define the eigenvalues of the S- and A-modes as cr^ = 
and <t^ = + iuj( A \ The spectral resolution of our calculations was 16 radial and 25 
latitudinal basis functions. The results were qualitatively confirmed by a number of runs 
with larger number of the basis functions. 



3. Results 

We calculate the dynamo solutions for a low level of the background turbulent 
diffusivity, choosing C v = 0.1, in order to approximately match the period of the eigenmodes 
with the solar cycle period. This corresponds to the background diffusivity ~ 10 8 m 2 s -1 . 
Figure [2] (left column) shows the linear-stability diagrams for the dynamo models with 
the meridional circulation speed values equal to Uq = 8, 12 and 16 ms -1 . The growth 
rate X^ = Re (cr^ A ^) of the first, most unstable dynamo mode (A-mode) is shown by 
the color-scale plots in the (C a , C&) plane, where C a and C$ are the free parameters that 
control the strength of the a and Ox J effects. We find that the dynamo instability 
region (represented by red color) changes significantly with the increase of the meridional 
circulation speed. For the slow meridional circulation it is found that the first A-mode is 
stable and steady in the absence of the a-effect (C a = 0). It has the excitation threshold of 
C$ « 0.013. In the opposite limit, when C$ = 0, the first mode is stable and oscillating. Its 




k 



Fig. 2. — Left column shows the linear-stability diagrams for the dynamo models with the 
meridional circulation speed Uq = 8, 12 and 16 m • s _1 from the top panel to bottom. The 

color scale shows the growth rate of the most unstable A-mode, \( A \ The contours show 

(o) 

77 

the oscillation frequency, uj^ a \ of this mode in the units of -55-. The right column shows 

R r . 



the relative growth rates of A- and S-modes: 



AA |A^)|-|A(*)| 



where S-mode has the 



A |A( A )| + |A(*)| 

symmetric toroidal field relative to the equator. The A-mode dominates in the regions 



colored in red. 



- 11 - 



oscillation frequency grows with the increase of the a-effect parameter C a . The excitation 

threshold is C a « 0.025 . The oscillation frequency of the mode at the threshold is about 

(o) 

R% 

The right column in Figure [2] shows the growth rate of the first A-mode relative to the 

|A^)|-U (S) | 

first S-mode. The relative difference is characterized by parameter w ', . This 

+ \\( s )\ 

helps to identify the regions in the parameter space (C a , Cs) where the A-mode dominates 
the S-mode. We find that in the case of Uo — 8 ms -1 the A-mode is dominant for C a > C$. 
In this regime we look for a solar-type dynamo solution, because in the solar dynamo the 
toroidal magnetic field is A- type (antisymmetric relative to the equator). 



For an example, we examine the case of Cs = 0, C a w 0.025 when the first A-mode 

(o) 

shows the snapshots of the magnetic field variation 



has the frequency « 8—^-. Figure 

it© 

inside the convection zone (top) and the time-latitude "butterfly" diagram for this mode 
(bottom). The snapshots show that the toroidal magnetic field is concentrated at the 
bottom of the convection zone, and the poloidal field is concentrated in the polar regions. 
Also, the toroidal magnetic field is globally distributed in the bulk of the convection zone. 
The maximum of the toroidal field distribution drifts to the equator at the bottom of the 
convection zone and moves to the pole near the surface. The bottom panel of Figure [3] 
shows the butterfly diagrams of the toroidal field at the bottom of the convection zone 
(color background) and for the radial magnetic fields at the surface (contour lines). The 
toroidal magnetic field evolution pattern has the polar and equatorial branches. The 
equatorial branch is strongly concentrated to equator. The phase relation between the 
radial magnetic field in the polar regions and the toroidal field in the equatorial regions 
is in agreement with observations of the polar magnetic field and the sunspot butterfly 
diagram, assuming that sunspots are formed by emerging toroidal magnetic field. However, 
this dynamo mode lacks the equatorial branch of the large-scale radial magnetic field, which 
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is also found in observations. The period of the dynamo period is about 12 years. This is 

as twice as short compared to the solar magnetic cycle. The period can be increased by 

further decreasing the diffusivity parameter C v by a factor of 2 (C^ ~ 0.05). However, this 

leads to decrease of the excitation threshold and increase of the effective magnetic Reynolds 
U R 

number (R M = — ). This means that the S-mode becomes dominant, and solution no 

Vt 

longer corresponds to the solar dynamo. Figure [2] shows that for the case of Uq = 16 ms _1 
the S-mode dominates. 

Inspecting the results in Figure 3 we can conclude that the increase of the meridional 
circulation speed has two main effects on the dynamo instability. Firstly, the larger Uq 
the larger the unstable area in the (C a , Cs) space, occupied by the first unstable A-mode 
(associated with a non-oscillating dynamo solution). Secondly, the S-mode becomes 
dominant near the instability threshold everywhere, both for the case of C§ = and 
arbitrary C a and for the case of arbitrary Cs and C a — 0. The combination of the a and 
Q x J effects can make the A-mode dominant, but it represents a non-oscillating dynamo 
solution. 

In another example, we examine the model, when the poloidal field is generated both by 

the a-effect and the Q x J-effect, e.g., C a = Cs = 0.015 and Uq — 8 ms" 1 . The oscillation 

(o) 

frequency of the first unstable A-mode is about 4— y. Near the excitation threshold, the 

R Q 

A-mode is highly dominant over the first S-mode. Figure [4] shows the snapshots of the 
magnetic field variations inside the convection zone (top) and the butterfly diagram for this 
mode (bottom). The snapshots of the magnetic field evolution inside the convection zone are 
similar to the previous case. However, the toroidal magnetic field is stronger concentrated at 
the bottom of convection zone, and the polar branch of the toroidal magnetic field evolution 
near is weaker near the surface. The period of dynamo wave is about 24 years, close to the 
solar cycle. Generally, we see a significantly better agreement with the observations than in 
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the previous case. 

In Figure [5] the dependence of the dynamo wave period on the speed of the meridional 



flow along the stability threshold shown. Contrary to previous results, e.g., (Bonanno et al. 



2002; Seehafer & Pipin, 2009) the period is not a monotonic function of the flow velocity. 
The main reason is that here we use the meridional circulation with a different depth 
dependence. We made a check how the presented results may depend on the distribution of 
the alpha effect. For this,, we switched off the effects of the turbulent mixing stratification, 



j^(u) = o in Eq.(10). The dynamo period as a function of the meridional flow speed for this 
case is shown in Figure fright). We find that dependence of the dynamo period on the 
flow speed is much stronger in case = for both types of dynamo. 



4. Discussion and conclusions 

We have studied kinematic axisymmetric mean-field dynamo models for a meridional 
circulation pattern with a deep seated stagnation point. This kind of circulation is suggested 
by the mean-field models of the angular momentum balance in the solar convection zone. 
We show that by adjusting the turbulent sources of the poloidal magnetic field generation 
and the turbulent diffusion strength it is possible to construct a mean-field dynamo model 
that resembles in some aspects the solar magnetic cycle. The most important features of 
the investigated models are the following. 

The maximum strength of the toroidal magnetic field, which is believed to be 
responsible for the sunspots production, is concentrated near the bottom of the convection 
zone. This field is transported to the equatorial regions by the meridional flow. The 
meridional component of the poloidal field is also concentrated at the bottom of the 
convection zone. The large-scale radial field is concentrated near the poles. It reverses sign 
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5 10 15 20 25 30 



TIME, [YR1 

Fig. 3. — The snapshots for the magnetic field variation inside convection zone (top) and 
the butterfly diagram for this mode(bottom), for the meridional flow speed Uq = 8 ms _1 and 
C a = 0.025 and C s = 0. 




5 10 15 20 25 30 35 40 



TIME, [YR] 



Fig. 4. — The as in Figure [3] for a dynamo with C a = Cs = 0.015 
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when the maximum of the toroidal field gets close to equator. This is not quite consistent 
with the solar observations, which show that the polar field reversals happen earlier in the 
cycle. Similar result was demonstrated in the kinematic flux-transport model by |Rempel 



(2006). His model had a qualitatively similar meridional circulation pattern, and the speed 
at the bottom of the convection zone was as twice as slow compared to our case. We 
believe that this feature (the phase relation) is inherent for this type of the meridional 



flow, which produced conveyor-belt like circulation of magnetic field (Dikpati et al. , 2004). 
The equatorward and poleward conveyor bands are not well connected in our case because 
circulation is quite weak in the bulk of convection zone. 

We find that including the combined action of the a and ft x J effects for the poloidal 
magnetic field generation improve the agreement of the model with observation. Contrary 
to the usual expectations that come from the results of the flux-transport dynamo model we 
find that the period of the dynamo cycle does not always become shorter when the speed 
of the meridional circulation increases. In our model this rule works for the amplitude of 
flow > 3ms -1 in the case of the a 2 Q dynamo with the a effect dependent on the density 
stratification, and for > 8ms -1 in the case of the a 2 SQ with the a effect dependent on both 
the density and the turbulent diffusivity stratifications. The dependence of the dynamo 
period on the flow amplitude is much stronger if the alpha effect does not depend on the 
turbulence intensity stratification, = 0. 

Thus, by measuring the distributions of the magnetic activity and meridional circulation 
characteristics on the Sun and possibly other cool stars we may get indirect information 
about contribution of the fix J effect to the dynamo and the relative contributions to the 
a effect due to density and the turbulent diffusivity stratifications as well. 

We conclude that the meridional flow pattern and speed have to be considered among 
the most important constrains on the stellar dynamo. Our results show the possibility of 



- 16 - 



using helioseismic observations of the meridional circulation for the diagnostic purpose of 
the solar dynamo, because the dynamo properties significantly depends on the depth of the 
flow stagnation point. 



5. Acknowledgements 



This work was supported by NASA LWS NNX09AJ85G grant and partially by RFBR 
grant 10-02-00148-a. 



- 17- 




6 8 

U , [m s" 1 ] 



Fig. 5. — Dependence of the dynamo period on the meridional flow speed Uq along the 
stability boundary of the most unstable dipolar mode for the dynamo models with the fix J 
effect (solid line) and without it (dashed line). Left, the results for complete a-effect and 



right for the a-effect with A w = (see, Eq.(lO)) 
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6. Appendix 

Here we describe the components of the mean-electromotive force which is used in the 
model. The tensor c^j represents the turbulent alpha effect, and in accordance with P08 it 
is given by 

= %{3^(A ( a) (e-A("))+/ 1 ( f(e.AW))}+ (10) 
+ ei e i {3^(/W(e.AW)+/W(e.AW))} 

3r?T { {e,kf + ei A?>) / 6 (a) + (e.kf + e.A^) / 8 (a) } , 



+ 

tensor 7^ describes the turbulent pumping 

Hi = {/i a) Ai p) + f[ a) (e • A^) e n } £m , - (11) 
and the 77^ term describes the anisotropic diffusion due to the Coriolis force and the x J 



effect (Radler, 1969), 



3 VT { (2/{ a) - e ijk - 2f[ a) ei e n e njk + /i%e fc } , (12) 



functions f^_ n y (given below) depend on the Coriolis number Q* = 2t c Qq and the typical 
convective turnover time in the mixing-length approximation is r c = l/u f . The turbulent 
diffusivity is parametrized in the form, t\t = C v r]!^\ where rj^ = is the characteristic 
mixing-length turbulent diffusivity, u f is the RMS convective velocity, £ is the mixing length, 
C v is a constant to control the intensity of turbulent mixing. The background turbulence is 
a state of turbulent flows in the absence of the mean magnetic fields and global rotation. 



This manuscript was prepared with the A AS IAT^X macros v5.2. 
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The others quantities in Eqs. 12h are: p is the density stratification scale, 



A w = Vlog (vt^ * s ^ e sca l e °f turbulent diffusivity, e = ft/ \ Q\ is a unit vector along 
the axis of rotation. Equations ( To|ll|12 ) take into account the influence of the fluctuating 



small-scale magnetic fields, which can be present in the background turbulence (see 



discussions 


in Frisch et al. , 


1975; Moffatt, 


1978; 


Vainshtein & Kit chat inov[ 1 1983; 


Kleeorin 


et al., 


1996; 


Brandenburg & Subramanian, 


2005 


b 2 

). In our paper, the parameter e = =, 

p pu 2 



which measures the ratio between the magnetic and kinetic energies of fluctuations in the 
background turbulence, is assumed equal to 1. This corresponds to the energy equipartition. 



Below we give the functions of the Coriolis number defining the dependence of the 
turbulent transport generation and diffusivities on the angular velocity: 

1 



(a 



/i 



+ 



f. 



A a 
Ji 



A a 

J6 

A a 
Ji 



A® 
J 10 

A" 
Jll 



n 



40* 2 
1 

ifF 2 
1 

60* 3 
1 



x arctan 
(0* 2 + 3) — 3 



0* 



((s-l)O* 2 + £-3) 



arctan O* 
0* 



3-s , 



30 



*3 



3 (O* 4 + GsQ* 2 + 10s - 5) ard g^* _ (( 8e + 5 )fr 2 + 3Q£ _ ^ 
3 (O* 4 + 3eQ* 2 + 5(s - 1)) arct ^ Q * _ ((4s + 5)0* 2 + 15(s - 1)) ) , 



1 f 3 ((3s - 11) O* 2 + 5s - 21) _ (4 ( £ - 3) O* 2 + 15s - 63) ) , 



480* 3 
1 

120* 3 



0* 

3 ((3s + 1) O* 2 + 4s - 2) arC ^ - - (5 (e + 1) 0* 2 + 12e - 6) ) , 

arctan 0* 



30 
1 



^ (3 (O* 2 + 1) (0- + s - 1) - ((2s + 1) 0* 2 + 3s - 3)) 



60* 3 
1 

20 
1 

20 



3 (0* 2 + 1) (0* 2 + 2s - 1) arcta * Q * - ((4s + 1) 0* 2 + 6s - 3) ) . 



({e + 1) 0* 2 + 3s - ((2s + 1)0- + 3s) 
^((20- + 3)-3(0- + l)^M 



